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Abstract. We revisit the Trans-Neptunian Objects (TNOs) color controversy allegedly solved bv lTegler & Romanishii]|i2003h . 

We debate the statistical approach of the quoted work and discuss why it can not draw the claimed conclusions, and reanalyze 
their data sample with a more adequate statistical test. We find evidence for the existence of two color groups among the 
Centaurs. Therefore, mixing both centaurs and TNOs populations lead to the erroneous conclusion of a global bimodality, while 
there is no evidence for two color groups in the TNOs population alone. We use quasi-simultaneous visible color measurements 
published for 20 centaurs (corresponding to about half of the identifie d objects of this class), and conclude on the existence of 
two groups. With the surface evolution model of Delsanti at ajj i2003h we discuss how the existence of two groups of Centaurs 
may be compatible with a continuous TNOs color distribution. 
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1. Introduction 

A large population of small icy bodies ex ists beyond t he or- 
bit of Neptun e. Fi rst speculated b v lLeonard (1930 ); Edgewortl] 
ill 9431 1 19491) and iKuiDeiU 19511). their existence was observa- 
tionally confirmed bv 'jewi tt & LutJ ( 1 1993 1) only eleven years 
ago. Considered as remnants from the formation of the so- 
lar system, these Trans-Neptunian Objects (TNOs) constitute 
the Edgeworth-Kuiper Belt (EKB), therefore also frequently 
known as Kuiper Belt Objects (KBOs). Currently, more than 
700 of them have been det ected. A different class o f objects, 
the Centaurs, was found bv lKowal & Gehrelsl l ll977 *); to date, 
more than 40 of these objects are known. Orbiting mainly be- 
tween Jupiter and Neptune, a strict dynamical definition does 
not exist and frequently some objects are both classified as 
Centaurs and Scattered Disk Objects (SDOs). Centaurs are be- 
lieved to be ex-TNOs in a tra nsition phase between the EKB 
and th e Jupiter-family comets ( Fernandez ' 198(f; 'D uncan et alJ 
Il988l) . Their exact origin location in the EKB has not yet been 
identified. WhUe it is currently believed Centaurs origin ate 
from the Scattered Disk Objects ('Du ncan & Levison ''l997^ — 
a "fuzzy" subclass of TNOs with highly eccentric and inclined 
orbits — an origin in the Plutinos (TNOs in a 3:2 mean mo- 
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tion orbital resonance with Neptune) has also been hypothe- 
sized ( Yu & Tremaine 1999). 

Since the very beginning of the first photometric measure- 
ments of these objects their visible color distribution alwa ys 
been very controversial. iTegler & Romanish in" ( 1998, '200(/) - 
hereafter TR98 and TROO, respectively - reported the identifi- 
cation of two separated color groups (blue-gray, i.e., solar like, 
and very red) from samples of 16 and 37 objects respectively. 
Other groups state d a continuous color s preading such as th e 
precursor wor k s bv | l-.uu & JewittI lll996l) ; ICireen et al. I lll997h ; 
Barucci et al. ( 1999), and subsequent studies (see the re- 
view by Doressoundiram 2Q04. for a summary). Jewitt & LutJ 
(200l|) analyzed the statistical significance of this bimodality 
on TR98 and TROO data samples and also on their own dataset 
of 28 objects. They concluded that both results were formally 
consistent with derivation from a uniform color distribution, 
i.e. un imodality. With 21 objects, measured in BVRIJ bands, 
iBarucci et al. ( 2001) reported the identification of four groups 
with a quasi-continuous spreading. However this grouping 
is seen only with V - I and V - J color information and 
becomes marginal w i th B - V and V - R colors. Finally, 
ITegler & RomanishinI ( l2003l) - hereafter TR03 - claimed the 
resolution of the TNOs color controversy, now with an enlarged 
sample of 55 objects (50 objects with measured B - V and 
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Fig. 1. B-V vs. V -R plot of TR03's data used to conclude for 
bimodality. Centaurs are highlighted and their bimodal behav- 
ior is obvious, contrarily to TNOs taken alone. 

V -R colors, and 5 with V -R only). Furthermore, with Monte 
Carlo simulations they showed that with typical error bars of 
0.05 magnitude other teams would not see any clear "gap" be- 
tween their two detected groups, concluding that observational 
methodologies leading to smaller error bars are the basis for 
these color groups findings. 

Whereas it is physically difficult to understand the exis- 
tence of two groups of colors, several explanations have been 
proposed to describe the continuous color r ange of TNOs. First, 
TNOs might have real intrinsic differences. ICTomeslll2003h pro- 
posed a migration model explaining the present composition of 
the EKB as the mixin g of bodies formed in very different parts 
of the Solar System. ILuu & JewittI ( 1 19961) explored a surface 
evolution mechanism: a competition between reddening space 
weathering and bluishing collisional resurfacing, which causes 
the continuous color spreading. However this model predicts a 
surface color variation with the objects rotation (e.g. within a 
few h ours) that has never been observed to date djewitt & Luul 
I2OOII) : it also implies color-or bital parameters correlations 
very d ifferent from observations dThebault & D oressoundira 
I2OO3I). Nevertheless, modeling is still in progress: Gil-Huttoi 
ll2002l) . with a more detailed space weathering process, pro- 
posed a new version of this mechanism but of difficult test- 
ing. Delsan ti at al. (2003) revisited the model bv .Luu & Je witt 
(11994 testing an additional resurfacing process: cometary ac- 
tivity, which is compatible with a surface homogeneity. 



2. Two color groups' analysis 

In the first work on TNOs bimodality (TR98), the two color 
groups seemed very clear. However, statistics were performed 
on too few objects to be really significant. Later, with an en- 
larged data sample, the separation between the two presumed 
color groups was much more tiny, even though still apparent 
(TROO). We should nonetheless notice that in Tegler and 
Romanishin's analysis Centaurs were mixed with TNOs. The 
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Fig. 2. B-V vs. V-R plot of TR03's full data. The "red cluster" 
of TNOs with q > 40 AU is included. 
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Fig. 3. B-V vs. V-R plot of all TNOs from TR03 excluding 
the Centaurs. Without the Centaurs, the absence of bimodality 
is clear. 



irradiative and collisional environment of Centaurs should be 
very different of that of TNOs, hence conclusions based on 
the colors of both populations altogether should be taken with 
great care. There are few statistical tests for the detection of 
more than one mode (group) in a distribution. Tests based on 
bins, as used by TR03, should be regarded with prudence as 
they are too d ependent on the bin's choice. We will use o nly 
the Dip Test ("Hartigan & Hartiganl fl98i: LHartigan"l985'), a 
distribution-free test that computes the maximum difference 
("dip") between the empirical distribution function and the 
unimodal distribution function that minimizes the differences. 

If we first restrict ourselves to the 32 objects with known 
B - R colors and q < 40 AU of TR03 — a sample for which 
they found two groups — , the dip test confirms a bimodality 
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Fig. 4. B - V vs. V - R plot of our analyzed Centaurs' sample. 
The bimodality seems clear and the statistics come out highly 
significant. 

for B - R colors with a significance level (S L) of 99.7%. 
However, if we remove the Centaurs from this sample we have 
bimodality only at SL ~ 55% for the 22 objects left, while 
for the 10 Centaurs alone we see it at SL = 94.7%. Even if 
the obtained SL is too low and the number of objects is too 
small to make strong conclusions, these results suggest that the 
bimodality found is most probably dictated by the Centaurs 
colors and can not be stated as general for TNOs (see figure[0. 
It should be noted also that, on TR03, these two groups were 
found without the inclusion of the "red cluster" of TNOs with 
perihelia distances above 40 AU. These "red cluster" objects 
with q > 40 AU are also characterized by a low inclination 
orbit and they most probably consti tute a separate group 
jLevison & Stern 1 ^2001': Brown '2001') However, they appear 
as a red cluster while compared to the full color distribution of 
all other TNOs, that continuously range from neutral to red. If 
TNOs with q < 40 AU divide in two groups one neutral and 
another red, we can no longer consider the "red cluster" as a 
separate group (based on color differences) since it is appar- 
ently equal to the red one found for q < 40 AU. Consequently, 
TR03 are implicitly dealing, in their work, with not two but 
three groups of objects: neutral and red TNOs with q < 40 AU, 
and red TNOs with <7 > 40 AU. When all objects are included 
in the analysis, the bimodality disappears (SL ~ 36%). If 
Centaurs are removed the significance level goes down to 
~ 9% (see figures |2l and |3}. Thus, while it does not seem 
justifiable to ignore the "red cluster" in the statistics, including 
those does "destroy" the evidence for bimodality, even with 
the inclusion of Centaurs. We will therefore investigate the 
possibility that the conclusion for two color groups is in fact 
being misguided by the superposition of a bimodal Centaur 
distribution over a continuous spread of colors of TNOs with 
q < 40 AU and a separate group of exclusively red TNOs with 
q > 40 AU or, globally, an unimodality for TNOs colors with 
an excess of red objects. 



We therefore took all Centaurs with available quasi- 
simultaneous observations of B - V, V - R and B - R colors, 
building a data sample of 20 objects which represents almost 
50% of all known Centaurs. Eight objects were obtained 
unde r the "ESQ Large Prog ram on Centaur s and TNOs", six 
from [ Peixinho et alJ (l2003h and two from iBoehnhardt et alJ 
(120021) a nd anot her six drawn from t he "Meudon Multicolor 
Survey" dPores soundira met al.ll2002l) . Si x more were added 
from t he compiled MBOSS database of iHainaut & Delsantil 
(l2002t) . since their average color indexes were computed from 
literature using only simultaneous magnitudes (see table 
While TR03 observational methodology also monitored the 
absence of significant magnitude variations, all their Centaurs 
were also independently measured by the previously quoted 
works. The dip test on our 20 objects data sample reveals 
bimodality for B - R colors with SL = 99.5%. Furthermore, 
by testing B-V and V - R colors separately, we see that the 
bimodal behavior is dominated by the V - R color index, since 
there is evidence for two V - R color groups at S L = 97.7% 
while for B-V alone it is inexistent. All dip test results 
are summarized in table |2l We conclude that Centaurs' 
population is indeed composed by two separate B - R color 
groups with strong evidence for a dominating V-R bimodality. 

In the color evolution model by IPelsanti at al.l ll2003l) . 
TNOs are resurfaced by a competition of 1) reddening by ir- 
radiation; 2) bluishing by non-disruptive collisions between 
TNOs; 3) bluishing by cometary activity (i.e. uniform re- 
deposition of neutral-colored dust Ufted by gas outbursts). 
Simulations using known Centaurs orbits lead to very red ob- 
jects regardless of the size. In the spaces between giant plan- 
ets, reddening by irradiation is the dominant process while 
collisions are scarce; cometary activity outbursts (triggered by 
collisions) are also possible but with a very low probability. 
Observed blue Centaurs colors (like for 2060 Chiron) are com- 
patible wi th a surfac e being recovered by the bound coma de- 
tected by Meechetal. (1997). It has been suggested that some 
short period comets and Cent aurs should be the fragm ents of 
objects ejected from the EKB jFarinella & D avis"* 1996). Such 
fragmentation would expose some fresh, volatile rich interior 
that would be compatible with the e xistence of active and/or 
blue Centaurs. IPelsanti a t al.' ('2003") concluded that very red 
Centaurs are old, fully irradiated objects (with a thick, dark 
irradiated mantle that prevents from spontaneous - i.e. non 
collision-induced - cometary activity outbursts), while blue 
centaurs might be fresh fragments recovering from an ejec- 
tion from the EKB. Unfortunately, their model does not cur- 
rently predict the proportion of blue vs red Centaurs. However, 
following this formalism, one can expect it is more likely to 
observe red Centaurs. We see with the sample used that the 
two groups (red/blue) are equally distributed. Assuming that 
the "excess" of blue objects is not an observational bias, we 
will have to invoke an iriiectiori of new blue Centaurs from 
the EKB. IYu & Tremaind d 19991) suggested that evolution of 
Plutinos onto Neptune-crossing orbit may dominate the flux 
of Jupiter-family comets. Consequently, Centaurs' population 
may be mainly populated by ex-Plutinos. An excess of blue 
Plutinos within the same size range as Centaurs is reported by 
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Table 1. Centaurs' colors 





Object 


B-V 


V 


-R 


B-R 




2002 GO," 


1.03+0.07 


0.75 


±0.04 


1.80 ±0.07 




2002 DHs" 


0.66 + 0.07 


0.47 


±0.05 


1.05 ±0.07 




I999XX143'' 


1.28±0.11 


0.65 


±0.07 


1.86 ±0.07 




1994 TA*" 


1.26 + 0.14 


0.67 


±0.08 


1.93 ±0.16 


(63252) 


2001 BL41" 


0.82 ± 0.06 


0.41 


±0.04 


1.18 ±0.06 


(60558) 


2000 ECqs" 


0.85 ± 0.08 


0.47 


±0.05 


1.32 ±0.08 


(55576) 


2002 GBio" 


1.09 ±0.05 


0.69 


±0.03 


1.77 ±0.05 


(54598) 


2000 QC243'' 


0.67 ± 0.03 


0.44 


±0.03 


1.11 ±0.03 


(52975) 


Cyllarus'' 


1.17 ±0.05 


0.71 


±0.04 


1.88 ±0.05 


(52872) 


Okyrhoe'' 


0.89 ±0.11 


0.43 


±0.08 


1.32 ±0.11 


(49036) 


Peliotf 


0.77 ±0.10 


0.47 


±0.10 


1.25 ±0.10 


(44539) 


1999 0X3" 


1.16 ±0.06 


0.74 


±0.07 


1.89 ±0.06 


(33128) 


1998 BU48'' 


1.01 ±0.04 


0.64 


±0.03 


1.65 ±0.04 


(31824) 


Hiatus'' 


1.05 ±0.02 


0.64 


±0.02 


1.69 ±0.02 


(10370) 


Hylonome 


0.77 ± 0.08 


0.38 


±0.06 


1.15 ±0.08 


(10199) 


Chaiiklo'' 


0.80 ± 0.05 


0.48 


±0.03 


1.28 ± 0.06 


(8405) 


Asbolus'' 


0.75 ± 0.04 


0.51 


±0.07 


1.26 ±0.08 


(7066) 


Nessus'' 


1.09 ±0.04 


0.79 


±0.04 


1.88 ±0.06 


(5145) 


Pholus'' 


1.30 ±0.10 


0.79 


±0.03 


2.09 ±0.10 


(2060) 


Chiron'' 


0.68 ± 0.04 


0.36 


±0.03 


1.04 ±0.05 



Colors are 



from, 




Hainaut & Delsanti 
DoresspundirarrL et jiL ( 20Q2) 

Table 2. Dip Test Results 



respec tively, "_JPgixinho_e £al] 12^ 
('2002'), ' iBoehnhardt et alJ 



Data Sample 



Color Dip" 



TR03 

... Centaurs+TNOs: q < 40 AU 
... Centaurs 
...TNOs: q< 40 AU 
...Centaurs+allTNOs 
... all TNOs 



32 
10 
22 
50 



B-R 
B-R 
B-R 
B-R 



0.1073 99.7% 

0.1389 94.7% 

0.0736 ~ 55% 

0.0442 ~ 36% 



40 B-R 0.0411 



^ 9% 



This work 
... Centaurs 
... Centaurs 
... Centaurs 



20 B-R 0.1211 99.5% 
20 B-V 0.0750 - 53% 
20 V-R 0.1156 97.7% 



" Number of objects, Measure of minimum difference from uni- 
modality, Dip's significance level 



IPeixinho et al.l ( l2003l) as also similar color-color correlations 
for Centaurs and Plutinos in opposition to SDOs correlations. 
These results give plausibility to our hypothesis. Nevertheless, 
different origins for the Centaurs population, or part of it, as an 
explanation of the two color groups can not be discarded. 

3. Conclusions 

In this work, we revisited the TNOs color controversy with a 
different and necessary twofold approach. First we argue that 
the dip test we used is a more appropriate tool to investigate 
the uni- or bi-modal distribution of colors. And second we 
suggest that the Centaur population has not to be mixed with 
TNO population in statistical studies since both populations are 
experiencing very different degrees of surface processing. We 
combined new and published datasets to show that the Centaur 



color distribution is bimodal while the TNO color distribu- 
tion is definitively unimodal. As Centaurs are presumed es- 
capees from the EKB , this new scheme is still compatible with 
evolution processes jPelsanti at aljEoOSi) . Moreover, our re- 
sults support a Plutino origin for the Centaurs. However, given 
the uncertainties on the transfer mechanisms of Centaurs from 
the EKB, we cannot rule out the possibility of true composi- 
tional diversity of Centaurs as an explanation of their color di- 
chotomy. 
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